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Abstract. We study, both theoretically and experimentally, the dynamical polarizability 
a{(o) of Rb 2 molecules in the rovibrational ground state of a^E+. Taking all relevant excited 
molecular bound states into account, we compute the complex-valued polarizability a{(o) 
for wave numbers up to 20000 cm^^ Our calculations are compared to experimental results 
at 1064.5nm (~ 9400cm^*) as well as at 830.4nm (~ 12000cm^^). Here, we discuss the 
measurements at 1064.5 nm. The ultracold Rb 2 molecules are trapped in the lowest Bloch band 
of a 3D optical lattice. Their polarizability is determined by lattice modulation spectroscopy 
which measures the potential depth for a given light intensity. Moreover, we investigate the 
decay of molecules in the optical lattice, where lifetimes of more than 2 s are observed. In 
addition, the dynamical polarizability for the state is calculated. We provide simple 

analytical expressions that reproduce the numerical results for a{(o) for all vibrational levels 
of a^E+ as well as X*E+. Precise knowledge of the molecular polarizability is essential 
for designing experiments with ultracold molecules as lifetimes and lattice depths are key 
parameters. Specifically the wavelength at ~ 1064 nm is of interest, since here, ultrastable 
high power lasers are available. 
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1. Introduction 

Owing to the extraordinary control over the internal and external degrees of freedom, ultracold 
molecules trapped in an optical lattice represent a system with many prospects for studies in 
ultracold physics and chemistry [1,2], the realization of molecular condensates [3], precision 
measurements of fundamental constants [4, 5, 6, 7] and quantum computation [8, 9] and 
simulation [10]. In the recent years, several groups have realized the preparation of optically 
trapped vibrational ground state (v = 0) molecules in either the lowest lying singlet or triplet 
potential [11, 12, 13, 14, 15]. Experiments with these molecules, e.g. ultracold collisions, 
are typically carried out in optical lattices or optical dipole traps [16, 17, 18]. In these 
environments, precise knowledge of the dynamical polarizability of molecules is important 
for well controlled experiments. 

The Rb 2 molecule is one of the few ultracold molecular species currently available, with 
which benchmark experiments for nonpolar molecules can be carried out. Here, we investigate 
the dynamical polarizability a (to) of a Rb 2 triplet molecule in the lowest rovibrational level of 
a^E+. A similar analysis for Cs 2 regarding the electronic ground state was previously 
carried out by Vexiau et al. [19]. In addition to calculations of the frequency dependent 
dynamical polarizability a (to), we present measurements of the real part Re{a(ft))} at a 
wavelength of A = 1064.5 nm. The experiments are performed with molecules trapped in 
the lowest Bloch band of a cubic 3D optical lattice which consists of three standing light 
waves with polarizations orthogonal to each other. By carrying out modulation spectroscopy 
on one of the standing light waves, we map out the energy band-structure for various light 
intensities. From these measurements Re{a(a))} is determined. Our experimental findings at 
A = 1064.5nm and also those at 830.4nm [13] agree well with the calculations. In addition, 
we experimentally investigate the decay time of the deeply bound molecules in a 3D optical 
lattice at 1064.5 nm for various lattice depths. Here, lifetimes of more than 2 s are observed. 
Furthermore, we present numerical results for the dynamical polarizabilities of the rovibronic 
ground state, i.e., the lowest rovibrational level of the X^E+ potential. For convenient 
application of our results, we provide a simple analytical expression and the corresponding 
effective parameters which can be used to reproduce the dynamical polarizabilities for all 
vibrational levels of both, the lowest singlet as well as triplet state outside the resonant 
wavelength regions. 

This article is organized as follows. In section 2, we give a brief, general introduction 
to the dynamical polarizability a((o) of a homonuclear diatomic molecule. Sections 3 to 
6 describe the calculations and measurements related to the polarizability of ^^Rb 2 in the 
rovibrational ground state of a^E+, along with a comparison of our results to reference values 
from literature. In section 7 we present calculations of the polarizability for the rovibrational 
ground state of X^E+. Afterwards, section 8 provides a simple expression which parametrizes 
the polarizability for all vibrational levels of the lowest singlet as well as triplet state. Tables 
with the corresponding parameters can be found in the Supplemental Material [20]. 
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2. Interaction of a diatomic molecule with light 


When a nonpolar molecule is subject to a linearly polarized electric field E = eEQCOs{(Ot) 
with amplitude Eq and unit polarization vector e, a dipole moment p =a {(o)E is induced. In 
general, a (to) is a tensor (see, e.g., [21]). For the sake of simplicity, we restrict ourselves to 
Hund’s case (b) molecules in the lowest rotational level = 0 of the nuclei, for which only the 
scalar isotropic polarizability a((o) is relevant (see, e.g., [21, 13]). Here, N = L-\-R, where L 
denotes the total electronic orbital angular momentum and R is the mechanical rotation of the 
atomic pair. The complex dynamical polarizability characterizes the response of a molecule 
to the electric field expressed by photon scattering and the ac Stark shift of molecular levels. 
This shift is directly linked to the dipole interaction potential 

U = -{p-E)l2 = -iRe{a(m)}£o" (1) 

and therefore to the real part of the polarizability Re{a(ft))} (see, e.g., [22]). In Eq. (I), the 
angled brackets (...) indicate time averaging. We note, that the dipole potential is attractive, 
when the sign of Re{a(a))} is positive and repulsive otherwise. The imaginary part of the 
polarizability Im{a((o)} is related to the power Eabs absorbed by the oscillator from the 
driving field, since 

^abs = {p-E) = ^mIm{a((o)}Eo. (2) 


We calculate the dynamical polarizability a((o) following the method described in Ref. [19]. 
The generic expression of the polarizability for a diatomic molecule in a state |/) is 
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Here, the angled brackets refer to the spatial integration over all internal coordinates of the 
system. The summation covers all the accessible dipole transitions with frequency tO/y and 
transition electric dipole moment d from the initial state |/) to final states |/) with line width 


r/- 


3. Calculation of a(o)) for a^E+ molecules 

3.E Relevant transitions 

If the molecule is initially in a vibrational level Va of the a^E+ state, all rovibrational levels 
(including the continuum) of the electronic potentials with ^E+ and ^Hg symmetry need to 
be accounted for in Eq. (3). As = 0 in the initial a^E+ state, only transitions towards final 
levels with total angular momentum 7 = 0,1,2 must be considered. Here, J = S-\-N, where 
S is the total electronic spin. Therefore, when considering a diatomic molecule it is usual to 
define two contributions to the isotropic polarizability a\ the parallel polarizability ct j along 
the molecular axis Z, which is related to dz, and the perpendicular polarizability a^, which 
is related to dx = dy. In general, a[| involves ^E+ —)■ ^E+ transitions and is related to 
3E+ —)■ transitions. One can show that a = (a|| + 2a^)/3 (see, e.g.. Refs. [19, 21]). 
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The expression given by Eq. (3) deals only with the transitions involving the two valence 
electrons of Rb 2 . Following Ref. [23], the contribution to the polarizability of the two Rb+ 
cores, hereafter referred to as etc, niust be taken into account, and is added to the results of 
Eq. (3). More details about this quantity are discussed in Appendix A.l. 

The first step of the calculations is to collect a set of accurate molecular potential energy 
curves (PECs) and transition electric dipole moments (TEDMs). The a^E+ PEC is taken 
from the spectroscopic study of Ref. [24]. For the excited molecular states and the related 
TEDMs from the a^E+ state, we use the same data as Refs. [25, 26], which we report in the 
Supplemental Material [20] attached to the present paper, for convenience. The PECs are 
displayed in Fig. 1(a), while the TEDMs are drawn in Fig. 10 (see Appendix A.2). These data 
are obtained by the quantum chemistry approach described in details in Ref. [27]. Briefly 
the Rb 2 molecule is considered as two valence electrons moving in the field of the two ionic 
Rb+ cores, which are represented by a large effective core potential (ECP) including a core 
polarization potential (CPP) [28, 29]. A full configuration interaction (FCI) is then performed 
on the two valence electrons, using a large gaussian basis set [30], with the CIPSI quantum 
chemistry code developed at Universite Paul Sabatier in Toulouse. It is worth mentioning 
that partial spectroscopic information is available on the 1^E+ state [31], the l^TIg state [32], 
and on the 2^11^ state [24], but no complete PEC has been extracted in these studies. As 
discussed for instance in Refs. [25, 32], the computed PECs are suitable to reproduce the 
observed data provided that they are slightly shifted in frequency (in terms of (0/{2tzc), by 
at most 100cm ^). We will estimate in section 5 the limited influence of such shifts on the 
results reported in the present work. Finally, the vibrational wave functions of levels |/) and 
I/) for the summation are obtained using the Mapped Fourier Grid Hamiltonian representation 
[33, 34]. 

3.2. Results 

The real and imaginary parts of the dynamical polarizability ay^=o{co) of a molecule in the 
vibrational ground state of the a^E+ potential are displayed in Fig. 1(b) and (c) as functions of 
the trapping laser frequency. The polarizabilities are expressed in atomic units (a.u.), which 
can be converted into SI units according to I a.u. = 4;reoaQ = 1.649 x lO^^^Jm^V^^, where 
uq denotes the Bohr radius and Cq is the vacuum permittivity. Note, for some applications, 
e.g., considerations related to the ac Stark shift, units of HzW^^cm^ (1 a.u. corresponds to 
4.6883572 X lO^HzW ^cm^) are advantageous. The sum in Eq. (3) has been truncated to 
include only the vibrational levels of the four lowest ^E+ states and the three lowest ^ITg 
states. Furthermore, electric-dipole-forbidden transitions are not considered in the sum, as 
they would appear as very weak and narrow resonances in the polarizability. The associated 
molecular data are collected in the Supplemental Material [20]. For simplicity, the natural 
lifetime Ty = iyf)^^ has been fixed to 10ns (jf ^ 2;r x 15 MHz) for all the excited molecular 
levels. 

Strongly oscillating patterns in both Re{avfl=o(w)} and Im{av^=o(®)} [see Figs. 1(b) 
and (c), respectively] correspond to frequency ranges of strong absorption which should 
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interat. dist. (ao) Re {ava=o} (a-u.) Im {ay^^o} (a-u.) 

Figure 1. (a) (solid blue lines) and (dashed red lines) potential curves of Rb 2 [25, 26]. 
The a^2i+ potential is drawn in black. In (b) and (c) the real and imaginary parts of the 
dynamical polarizability oc,.^=o of Rb 2 molecules in the rovibrational ground level (va = 0, 
R = 0) of the a^E+ molecular state are shown as a function of (o/(2nc). The two wavelengths 
used in the experiments are indicated by dashed horizontal lines. Furthermore, the red dashed 
vertical line in (b) represents zero polarizability. 


be disregarded for trapping purpose. The real part smoothly inereases from the statie 
polarizability av^=o{co = 0) = 698.5 a.u. up to the bottom of the l^S^ potential well, reaehing 
3147 a.u. at the wavelength of the trapping laser used in the present experiment (1064.5 nm). 
In the same region the imaginary part increases from about 10 ^ a.u. at CO = 0 to 10~^ a.u. at 
the trapping laser frequency which leads to a correspondingly larger photon scattering rate. 

It is difficult to provide a well-defined error bar on the theoretical values of the dynamical 
polarizabilities as their accuracy depends on the considered wavelength. Various causes of 
global inaccuracies have been analyzed in depth in Ref. [35]. First, the choice of a constant 
radiative lifetime for all excited levels in Eq. (3) influences only the strongly oscillating 
regions of the polarizability, changing the amplitude of the resonances. We checked that this 
approximation has no effect in the smoothly varying regions which are relevant for trapping 
experiments. We verified also that adding a couple of upper electronic states in the sum of 
Eq. (3) contributes to the polarizability for less than 1%. Moreover, the first excited E and IT 
states contribute together for more than 90% to the polarizability. Usually, these are the most 
well known states either because accurate spectroscopic results are available, or because they 
are well-determined by quantum chemistry calculations. The accuracy of TEDMs is tedious 
to analyze as their experimental determination relies on line intensities which are difficult to 
measure accurately. However, one argument in favor of the accuracy of the TEDMs results 
when comparing the values obtained from different methods. Eor instance, in Refs. [36, 37], 
with respect to various alkali-metal dimers, the TEDMs computed by two different methods 
are found to agree within 2%. Einally, an indication for the accuracy of the present work 
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is provided by the measurement of the dynamical polarizability for the v = 0, 7 = 0 level 
of the Cs 2 electronic ground state at 1064.5 nm, which is quite far away from the lowest 
resonant region [19]. The experimentally determined value with respect to A = 1064.5 nm 
is 2.42(15) X Ctcs [39], where OCqs is the dynamical polarizability of the Cs atom. This is in 
remarkable agreement with the computed value of 2.48 x ocqs [19]. 

4. Measurement of Re{a(fD)} for 

4.1. Experimental setup and measurement scheme 

The experiments presented in this work are carried out with a pure sample of about 1.5 x 10^ 
^^Rb 2 molecules prepared in the rovibrational ground state of the a^E+ potential and trapped 
in a 3D optical lattice. There is no more than a single molecule per lattice site and the 
temperature of the sample is about 1 /iK. As described in detail in Refs. [13, 21, 38], the 
molecules are prepared as follows. An ultracold thermal cloud of spin-polarized *^Rb atoms 
(/Rb = 1’ Rb ^ adiabatically loaded into the lowest Bloch band of a 3D optical lattice 
at a wavelength of A = 1064.5 nm. The lattice is formed by a superposition of three linearly 
polarized standing light waves with polarizations orthogonal to each other, see Fig. 2(a). The 
three lattice beams are derived from the same laser source with a linewidth of a few kHz and 
have relative intensity fluctuations of less than 10^^. In order to avoid interference effects, 
the frequencies of the standing waves are offset by about 100 MHz relative to each other. 
At the location of the atomic sample the beam waists (1/e^ radii) are about 130/im and 
the maximum available power per beam is about 3.5 W. By slowly crossing the magnetic 
Feshbach resonance at 1007.4 G we produce weakly bound diatomic molecules. After a 
purification step which removes remaining atoms, a STIRAP (stimulated Raman adiabatic 
passage) is performed at 1000 G, transferring the dimers into the rovibrational ground state 
(Vfl = 0, N = 0, m/v = 0) of the a^E+ potential. Here, m^v is the projection of N on the 
quantization axis defined by the direction of the magnetic field B [cf. Fig. 2(a)]. The 
molecule has positive total parity, total electronic spin 5=1, total nuclear spin 1 = 3 and 
is further characterized by the quantum number f = 2 (f = 5 + 1). Moreover, the total angular 
momentum is given hy F = 2 (F = f + J) and its projection is mp = 2. Henceforth, we simply 
refer to these molecules as ”va = 0 molecules”. 

According to Eq. (1), Re{a(a))} = A\U\/Eq, i.e., the real part of the dynamical 
polarizability can be determined by measurements of the potential depth \U\ and the electric 
field amplitude Eq of an optical trap. For the case of a cubic 3D optical lattice with orthogonal 
polarizations, the trapping potential is given by V {x,y,z) = L/3=xj,z^j3(jS) where the 

V^m = -\V\fCos^(k^ + i^ls) ( 4 ) 

represent the contributions of the standing waves of directions /3 = x^y^z with k = 2%/X 
being the wave number of the lattice beams. We now only consider the part of the lattice in 
the vertical z direction since this axis is the only one relevant for the measurements of the 
potential depth in the present work. Therefore, we define ^^(z) = y(z), \U\z= \U\ and 0z = ^. 
The phase 0 is a function of the laser wavelength A because the standing light wave is created 
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Figure 2. (a) Experimental scheme. Three retroreflected laser beams with polarizations 
orthogonal to each other form a 3D optical lattice. A trapped cloud of Rb 2 molecules is 
sketched in the intersection of the laser beams. The optical lattice is located between two 
Helmholtz coils which create a magnetic held B. In our setup, the direction of B represents 
the quantization axis, (b) Illustration of amplitude (left) and phase (right) modulation 
spectroscopy. 


by retroreflecting the laser beam from a fixed mirror at position Zm- It is given by ^ = Anzm!^ 
at z = 0. 


4.2. Lattice modulation spectroscopy 

In order to obtain |t/|, we earry out lattice modulation spectroscopy (see, e.g., [12,40,41,42]). 
For this, we either modulate |t/| by periodically changing the intensity of the standing light 
wave (amplitude modulation) or we modulate the phase <j) by periodically changing the laser 
wavelength A (phase modulation). Resonant amplitude (phase) modulation drives transitions 
from the lowest Bloch band (n = 0), in which the molecules have been initially prepared, 
to even (odd)-numbered excited lattice bands [see Fig. 2(b)]. This can cause either direct 
loss from the trap owing to heating (see, e.g., [40]) or molecules in higher lattice bands 
collide with each other and those of the lowest Bloch band, respectively, resulting in decay 
to nonobservable states. In consequence, resonant excitation leads to a decreased molecular 
signal in our measurements. 

For amplitude modulation spectroscopy, we modulate the intensity of the lattice laser 
beam sinusoidally by a few percent. When performing phase modulation spectroscopy, we 
modulate the laser frequency by a few MHz corresponding to a phase difference on the order 
of a few lO^^rad as Zm ~ 0.4 m. The modulation duration is typically on the order of 1 ms. 
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frequency v (kHz) 

Figure 3. Amplitude [(a) and (b)] and phase (c) modulation spectra of weakly bound Feshbach 
molecules (a) and molecules in the rovibrational ground state of the a^E+ potential [(b) and 
(c)]. We measure the fraction of remaining molecules N/Nq as a function of the modulation 
frequency v. The statistical error of each data point is in the range of ±(0.05 — 0.15). Here, 
the numbers Nq are given by the asymptotic limits of Lorentzian fits (solid blue lines). The 
resulting center frequencies of the resonances are illustrated as black vertical lines, while the 
red dashed lines indicate the corresponding uncertainties. We note, that the mean intensity of 
the lattice beam used for modulation in (b) is 30% less than in (a) and (c). 


At the end of each experimental cycle (which takes about 40 s), the remaining number N 
of molecules is measured. We only find molecules in the lowest Bloch band, not in higher 
bands. In order to determine the molecule number, we reverse the STIRAP and dissociate the 
resulting Feshbach dimers by sweeping over the Feshbach resonance. Then, the generated 
atoms are detected via absorption imaging. By comparing the resonant transition frequencies 
observed in the modulation spectra to the energy band-structure of the sinusoidal lattice, the 
lattice depth \U\ is deduced. This will be explained in detail further below. 

Figure 3 shows measured excitation spectra of Feshbach (a) and Vq = 0 (b,c) molecules, 
obtained via amplitude or phase modulation spectroscopy. A single data point typically 
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consists of 5 to 30 repetitions of the experiment (For a given speetrum the number of 
repetitions is eonstant). Fig. 3(a) as well as (b) exhibit a prominent resonanee after amplitude 
modulation. This resonanee is related to a transition from the lowest Bloeh band (n = 0) to the 
seeond exeited lattiee band (n = 2). Speetrum (a) for Feshbaeh moleeules in addition shows 
a broad shoulder at around 50 kHz whieh we attribute to a resonant transition from n = 0 to 
n = 4. Due to the large width of this resonanee the uneertainty in the determination of its 
eenter frequeney is relatively large. Speetrum (b) in Fig. 3 also features a seeond resonanee 
dip, but here it is loeated at about half the frequeney of the prominent one. This resonanee 
dip ean be assigned to a transition from the lowest Bloeh band to the seeond exeited band, 
involving two identieal “quanta” with frequeney V. It is known (see, e.g., [40]) that sueh 
subharmonie resonanees exist. To be eonsistent, a similar sub-harmonie resonanee dip should 
be present in Fig. 3(a) at about 15 kHz (as indieated by the vertieal, dashed arrow). Indeed, at 
that position the data points seem to be systematieally below the fit eurve with respeet to the 
prominent peak. However, the eorresponding signal (if at all) is very weak, partially due to its 
position at the steep flank of the prominent resonanee. 

Now, we turn to Fig. 3(o), whieh shows an exeitation speetrum after phase modulation for 
Va = 0 moleeules. We observe two resonanees of similar strength, both of whieh we assign 
to the transition from n = 0 to n = 1. The dip at lower frequeney is again a subharmonie 
resonanee. Surprisingly, it is stronger than the harmonie one at about 25 kHz. We attribute this 
to a purely teehnieal issue, as the strength of the phase modulation varied with the frequeney 
in our setup. However, we have verified the assignment of the resonanees by eomparison to 
the eorresponding amplitude modulation speetra. 

We ealeulate the Bloeh bands by diagonalizing the Hamilton operator for the lattiee in 
ID (negleeting gravitation). 


H = - 


<92 


+ ^(z), 


(5) 


2m dz^ 

whieh is partieularly simple in momentum spaee (see, e.g., [43]). Here, m is the mass of 
a moleeule, i.e., twiee the mass of a atom. Figure 4 shows the ealeulated energy 
eigenvalues as a funetion of the lattiee depth. The energies are given in terms of the 
reeoil energy Er = h^/{2mX^), with h being Planek’s eonstant. As we do not speeify the 
quasimomentum, the energy eigenvalues form bands whieh are broad for low lattiee depths. 
However, the bands n = 0ton = 2 are quite narrow for lattiee depths above ~ 40Er. This is the 
regime where we take most of our measurements. Having measured the resonant exeitation 
frequeneies after modulation we eould in prineiple use Fig. 4 to read off the eorresponding 
lattiee depth \U\. We refine this method and at the same time eheek for eonsisteney as follows. 

In the experiment we eontrol the lattiee depth \U\ via the laser beam power P that ean 
be measured using photodiodes. The square to the eleetrieal field Eq is proportional to P. 
Consequently \U\<^ P, i.e., the preeise value of |t/| is known up to a ealibration faetor (whieh 
depends linearly on the dynamieal polarizability). Thus, given a moleeular state, we should 
be able to adjust the ealibration faetor sueh that all data obtained for various powers P mateh 
the band strueture ealeulation. The measured data points in Fig. 4 elearly show that this works 
quite well, both for deeply bound moleeules (red) and Feshbaeh moleeules (blue). In this 
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Figure 4. Energy band-structure with zero energy corresponding to the center of the 
lowest Bloch band. Solid lines are calculations for a single lattice direction, where the 
numbers 1. to 8. give the band index n. The data points are obtained excitation energies 
for X = 1064.5 nm, stemming from amplitude (circles) or phase (squares) modulation 
spectroscopy. Red (blue) plot symbols indicate measurements for triplet rovibrational ground 
state (Feshbach) molecules. Here, the experimental results are shown after independently 
htting the data for each molecular species to the band-structure calculation. By doing so, we 
determine the individual calibration factors and therefore the lattice depths \U\ (see also text). 
The horizontal etTor bars represent the resulting uncertainties of |t/|, whereas the vertical error 
bars are given by the uncertainties of the Lorentzian hts in the modulation spectra. 


procedure we do not aeeount for the transitions from n = 0 to n = 4 owing to the large 
uneertainties of the eorresponding resonanees in the exeitation speetra. Nevertheless, these 
data points are shown in the plot for eomparison. 

In addition to \U\, the eleetrieal field amplitude Eq of the optieal lattiee has to be 
determined in order to infer the dynamieal polarizability «((») [see Eq. (1)]. We ean 
eireumvent this by refereneing the measurements on the lattiee depth \U\ for the moleeules in 
the rovibrational ground state of the a^E+ potential to similar measurements with Feshbaeh 
moleeules, of whieh the polarizability aFesh(w) is known to be twiee the one of a Rb atom 
aRb(tt)) in the electronic ground state [19]. Aeeording to Eq. (1) the lattiee depths |t/va=o| and 
|f^Fesh| for the Vfl = 0 and Eeshbaeh moleeules are related by 

\^Va=0\ _ Ro(Ogva=0) 

|f^Fesh| Ro(OCFesh) 

for a given lattiee beam intensity, i.e., a given Eq. 

Erom our experiments at A = 1064.5nm we obtain Re {ava=o} = (2.5 ±0.1) x Rejapesh}. 
whereas Re{avfl=o} = (0.1 ±0.02) x Rejapesh} was found at A = 830.4nm [13]. Using our 
ealeulated atomie values of 685.8 a.u (1064.5 nm) and 2995.9a.u (830.4nm) yields moleeular 
polarizabilities Re{avfl=o} of 3430 ± 140 a.u (1064.5 nm) and 600 ± 120 a.u (830.4 nm), 
respeetively. 
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5. Comparison of results 

Table 1 shows our measured and ealeulated polarizabilities along with results of other 
referenees. First, it should be noted that our theoretieal atomie polarizabilities (ineluding 
only the 5s —5p and 55 — 6p transition frequeneies from the NIST database [44] at 828.4 nm 
and at 1060.1 nm are in good agreement with the ones of Refs. [45, 46] whieh eonsider the 
5s —5p transition frequeney from the NIST database, and ab initio values for the frequeneies 
up to the 55 — Sp transitions. 

We find good agreement between our theoretieal (3147 a.u.) and experimental (3430 ± 
140 a.u.) results for the moleeular polarizability Re{ava=o} at 1064.5 nm (see table 1). In 
eontrast, the agreement for the polarizability at 830.4 nm of our former measurements [13] 
(600 ± 120 a.u.) with the present ealeulations (875.8 a.u.) is somewhat poor. In view of 
this diserepaney we want to estimate the influenee of slight shifts of the potential energy 
eurves on the ealeulations. The potential well depths of the 1^S+ and the l^TIg states used 
in the eomputation are smaller by 72 em^^ and 51 em^^ with respeet to the experimental 
determinations of Refs. [31] and [32], respeetively. Sueh shifts would lead to a ehange 
in the ealeulated polarizability of about 10% for the partieular wavelengths of 830.4 nm 
and 1064.5 nm. This sets a range for the uneertainty of the ealeulated polarizability that 
arises from the uneertainty of the PECs. Furthermore, in terms of the experiments, we 
note that in Ref. [13] a different method to determine the dynamieal polarizability was used. 
The polarizability was inferred from the oseillating dynamies of moleeular wave paekets 
that oeeurred when v = 0 moleeules were suddenly loaded into several Bloeh bands of the 
optieal lattiee. This leaves potentially room for a systematie diserepaney between the two 
measurements. With respeet to the statie polarizability, i.e., CO = 0, the ealeulations presented 
in this work give 698.5 a.u. (ef. table 1) for the = 0 moleeules. This value aetually agrees 
well with the one previously reported in Ref. [51], 677.5 a.u., as it was not ineluding the 
eontribution of etc = 2 x a(Rb+) = 18.2a.u. [47]. 

Figure 5 is a zoom into Fig. 1(b) showing the ealeulated real part of the dynamieal 
polarizability of a = 0 moleeule (solid blaek lines). In addition, Re{a(co)} for a Feshbaeh 
moleeule is plotted (dashed red lines), whieh is given by twiee the atomie polarizability. 
The two wavelengths used in our experiments (830.4 nm and 1064.5 nm) are indieated as 
vertieal green dashed lines. Outside the resonant and therefore lossy regions in Fig. 5 (vertieal 
blaek bands) the two polarizability eurves never eross. Thus, there is no so-ealled ’’magie” 
wavelength, where the ae Stark shift of the two moleeular states eaused by the trapping light 
is equal. Sueh state-insensitive trapping eonditions ean be benefieial, e.g., when eonverting 
Feshbaeh moleeules to deeply bound states. Speeifieally, in Ref. [13], owing to the large 
differenee of the dynamieal polarizabilities at A = 830.4 nm, the STIRAP transfer of Rb 2 
from the Feshbaeh level to = 0 populated several lattiee bands. 

We have again studied this issue in this work and find that population of higher lattiee 
bands ean be suppressed even in the absenee of a magie wavelength when working with deep 
lattiees. At A = 1064.5 nm there is still a faetor of 2.5 differenee in polarizability between 
Feshbaeh and = 0 moleeules. For an initial (final) lattiee depth of 50 Er (125 Er) at 
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Table 1. Measured and calculated polarizabilities Re{a} in a.u. for atoms and ^^Rb 2 
molecules in the rovibrational ground state of a^Z+. Here, the abbreviations ”tw exp (theo)” 
mean ’’this work, experimental (theoretical)”. 


Species 

A (nm) 

Re{a} (a.u.) 

Ref. 

00 

828.4 

3132±3 

[45] 

00 

828.4 

3131.4 

tw theo 

00 

830.4 

2995.9 

tw theo 

*^Rb2 

830.4 

875.8: 

tw theo 

*^Rb2 

830.4 

600 ±120 

[13], using ctRb = 2995.9 a.u. 

00 

1060.1 

692.7 

tw theo 

00 

1060.1 

693.5 ±0.9 

[46] 

00 

1064.5 

685.8 

tw theo 

00 

fO 

1064.5 

3147 

tw theo 

^^Rb2 

1064.5 

3430 ±140 

tw exp, using ctRb = 685.8a.u. 

^^Rb2 

1064.5 

3200 ±500 

[48] 

00 

oo 

318.6±0.6 

[50] 

00 

oo 

317.9 

tw theo 

*^Rb2 

oo 

698.5 

tw theo 

^^Rb2 

oo 

677.5 

[51] 


1064.5 nm a calculation of the wavefunction overlap for the Bloch states shows that still 97% 
of the population stays in the lowest Bloch band after the STIRAP. In addition, we are able 
to energetically resolve the lattice bands during STIRAP as n = 0 and n = 2 are separated 
by about 40kHz at \U\ = 125 Er (cf. Fig. 4). This strongly increases the selectivity of the 
transition (see, e.g., [52]). Indeed, in our experiments we do not observe any significant 
population of higher bands. 

As can be seen in Fig. 5, the absolute dynamical polarizability of the triplet rovibrational 
ground state molecules at 1064.5 nm is about four times larger than at 830.4 nm. This 
is convenient since it results in a four times deeper interaction potential at the same 
laser intensity. For longer wavelengths than 1064.5 nm. Fig. 5 reveals, that the dynamical 
polarizabilities of = 0 molecules and Feshbach molecules approach each other. Hence 
working at even longer wavelengths than 1064.5 nm might be advantageous for some 
applications. 

6. Lifetime of the molecules 


According to Eq. (2), the imaginary part of the dynamical polarizability, Im{a}, is linked to 
the light power absorbed by a molecule, Pabs> which in turn can be expressed in terms of the 
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Figure 5. Real parts of the dynamical polarizabilities of triplet rovibrational ground state 
molecules (solid black lines) and Feshbach molecules (dashed red lines). The circles represent 
the experimental results given in table 1 with the corresponding wavelengths indicated by 
green dashed vertical lines. 


photon scattering rate Tgc = /i ^Pabs [22]. Using this and Eq. (1), Esc can be written as 

2U lm{aico)} 

h Re{a((0)}' ^ ^ 

For a 3D optical lattice with equal lattice depths \U\ in each direction the scattering rate 
is given by = 3rsc. As an example, we consider the case of \U\ = 50 Er at X. = 
1064.5 nm. Then, the corresponding values for the polarizability obtained in the present 
work, Im(a) = 0.96 x 10^^ a.u and Re{a} = 3430a.u., yield Esc = 0.18 s^^ Note, this 
calculation only accounts for an ideal optical lattice. As there is always background light that 
does not contribute to the lattice, the estimated value for the scattering rate represents just a 
lower bound. Once a photon is absorbed, the molecule is excited and typically decays to a 
nonobservable state. Assuming excitation to be the only loss-mechanism a lifetime T = 1.9 s 
of the Vfl = 0 molecules is expected in a 50 E/? deep 3D optical lattice at 1064.5 nm. 

We experimentally investigate the lifetimes of the molecules in the rovibrational ground 
state of a^E+ by varying the holding time in the lattice. Figure 6 shows lifetime 
measurements of Vq = 0 molecules for various potential depths \U\, which are adjusted to 
be equal in each direction. Applying an exponential fit, we obtain a 1/e decay time T of more 
than 2 s for both our measurements at \U\ =33Er and \U\ = 5SEr. 

In order to estimate possible loss induced by inelastic molecular collisions, we calculate 
the tunneling rates Etu between adjacent lattice sites within the lowest Bloch band (n = 0). 
When considering a lattice depth of 33 Er (5SEr) one obtains E^ = 2.37 s^^ (Etu = 0.09 s^^). 
In our setup at most 20% of the lattice sites are occupied in the region of highest molecule 
density. Thus, for \U\ = 33 Er decay due to collisions cannot be neglected, whereas for 
\U\ = 5SEr and beyond the only relevant loss mechanism is photon scattering. For such 
deep lattices, the lifetime T scales directly inversely with the lattice depth \U\. We confirm 
this for the measurements at |t/i| = 5SEr and |f72| = 125since the ratio of the lifetimes 
Ti/t 2 = 2.20 is close to |t/ 2 |/|t/i| = 2.16. 
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Figure 6. Decay of triplet rovibrational ground state molecules trapped in a 3D optical lattice 
at 1064.5 nm with equal potential depths \U\ in each direction. Shown is the fraction of 
remaining molecules N/Nq as a function of the holding time fj, in the lattice. Square plot 
symbols (red circles) correspond to a lattice depth of 125 Er (5^Er). Each data set typically 
consists of 10 to 15 repetitions of the experiment, where the statistical error of a data point is 
on the order of ±0.1. Solid lines are exponential hts to the data. Here, Nq is given by the values 
of these fits at fh = 0 s. In general, the absolute molecule numbers No are about 1.5 x lO'^. The 
inset depicts the resulting 1/e decay times for various lattice depths. 


7. Polarizability ofX^E+ molecules 

As the agreement between ealeulations and measurements for the triplet moleeules is in 
general good, we also provide ealeulations for the singlet ground state of ^^Rb 2 . Using the 
same approaeh as above (see also [19]), we eompute the dynamieal polarizability ay^=o{a}), 
i.e., with respeet to the vx = 0, J = N = 0 level. The X^E+ PEC has been derived from 
speetroseopie data of Ref. [53]. The A^E+ and the b^n„ PECs and the related spin-orbit 
eoupling between those two states are taken from Ref. [54]. The PECs for all the other states 
and for the TEDMs are taken from the eomputations reported in Refs. [25, 26]. Again, the 
sum in Eq. (3) has been truneated to inelude only the levels of the four lowest ^E+ states 
and the three lowest states. The natural lifetime of the exeited levels has been fixed at 
10ns. Results are presented in Eig. 7, showing that two magie wavelengths ean be identified 
at 990.1 nm and 1047.2 nm. The latter is loeated elose to the region of strong absorption 
resonanees and eonsequently, from the imaginary part [ef. Eig. 7(b)], the photon seattering 
rate at 1047.2 nm is expeeted to be about four times larger than at 990.1 nm. 

8. Parametrization of the polarizability 

8.1. Rovibrational ground state ofa^'L^ 

In general, using figures (e.g., Eigs. 1 and 7) to read off the dynamieal polarizabilities 
at speeifie wavelengths is eumbersome. Therefore, we provide here a simple analytieal 
fitfunetion and parameters that allow for reprodueing the numerieal results with respeet to 







15 




Figure 7. Real part (a) and imaginary part (b) of the dynamical polarizability of ^^Rb 2 
X*X;+ molecules in their vx = 0, = 0 level. In (a), for comparison, also the numerical 

results corresponding to 2Re{aRb} are shown (red dashed line), representing the real part of 
the dynamical polarizability of Feshbach molecules. Two magic wavelengths, where both 
polarizabilities are equal, are obtained and indicated by vertical dashed lines. 


nonresonant wavelength regimes. In the infrared and optieal domain we are studying, the 
main eontributions to the polarizability of the X (a) state outside of resonanees eome from the 
transitions towards the first exeited state and the first state. Thus we attempt to 

model the polarizability by redueing those transitions to a single effeetive transition towards 
eaeh of the two different symmetries. The approximate real part of the polarizability is 
expressed as 


«eff(®) 


2(0e 




-dy 


(0^ 


+ 


2(On 




(0^ 


:^n 


+ OCf 


( 8 ) 


with coz (resp. cun) the effeetive transition frequeneies and dz (resp. du) the eorresponding 
effeetive dipole moments. We have isolated in this expression the eore polarizability etc as its 
frequeney dependenee is mueh weaker than the one of the terms eoming from valenee eleetron 
exeitation (see Appendix A.l). The imaginary part of the polarizability is negleeted sinee the 
model is designed for the ranges outside the resonant regions. 

We extraet the effeetive parameters from a fit to the full numerieal results using Eq. (8). 
The results with respeet to the a^E+ and X^E+ rovibrational ground state are shown in 
Fig. 8. For the triplet ease, data points with frequeneies elose to resonanees, i.e., from 
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Ci}/{lKc) = 9527 to llOlScm^^ and above 13173 cm^^, are exeluded from the fit. We 
obtain (Oz/i'lKc) = 10112.33em^\ dz = 2.792881 a.u., (0n/(2;rc) = 13481.32om^^ and 
du = 3.211023 a.u. Note, in terms of dipole moments, atomie units ean be eonverted into SI 
units aeeording to 1 a.u. = cuq = 8.478 x 10^^® Ams. Using these parameters, the effeetive 
polarizability reproduees the numerieal results in the fitted region to within a relative root 
mean square value (rRMS) of around 1% [see Fig. 8(a)]. The rRMS value is defined by 


rRMS = W ) — 



«eff(tQ/) - a{(Oi) 


(9) 


with a{(Oi) being the numerieal values, and aeff(tOi) the fitted ones. Here, M is the number 
of eonsidered values ft),, whieh depends on the vibrational level as the resonant frequeney 
regions exluded from the fit vary. We point out that if we take the transition dipole moment dz 
(resp. dx = dy) at the equilibrium distanee of the a^E+ state and multiply it by the appropriate 
Honl-London faetor (1/ for E+ — E+ transition and a/2/3 for E+ — H transition), we get 
the value 2.767 a.u. (resp. 3.142a.u.), very elose to the effeetive dipole moment found above. 
Moreover, the effeetive transition frequeneies eorrespond roughly to the average frequeneies 
of transitions with favorable Franek-Condon faetor. 


8.2. Rovibrational ground state ofX^'L^ 

A similar fit ean be performed for the dynamieal polarizability of X^E+ moleeules in the (v^ = 
0, A = 0) level [Fig. 8(b)]. Frequeney domains from (o/{2tzc) = 9432om^^ to 9960em^^, 
from 10613om^^ to 12890om^^ and above 14736em^\ eorresponding toresonanees towards 
the b3n„,AiE+andBin„ exeited states, respeetively, were exeluded from the fit. We omitted 
to aeeount for levels of the b^n„ state as they have a very small singlet eharaeter. Then, the 
fit to the numerieally ealeulated dynamieal polarizability yields (Oy/ {2tzc) = 11450.31 em^^, 
= 2.647731 a.u., ft)n/(2;rc) = 15019.57 em^^ and Jn = 2.965774 a.u. Again, we point out 
that the effeetive transition dipole moments obtained are very elose to the eleetronie transition 
dipole moments taken at equilibrium distanee multiplied by the Honl-London faetors, i.e., 
2.613 a.u. and 2.959 a.u., respeetively. Here, the rRMS of the fit is 0.5%. 

8.3. Excited vibrational states 

Both, the frequeney domains with good Franek-Condon faetors and transition dipole moments 
depend signifieantly on the initial vibrational level. This is refieeted in the variation of the 
effeetive parameters when we perform an individual fit for eaeh vibrational level (with N = 0) 
of the a^E+ (va = 0 to 40) and X^E+ (yx = 0 to 124) states. The eorresponding parameters, 
whieh ean be used to reproduee the dynamieal polarizabilities in the nonresonant frequeney 
domains, are reported in the Supplemental Material [20]. Figure 9 shows the resulting rRMS 
values as a funetion of the vibrational level with respeet to the lowest triplet and singlet state. 

Using the ansatz of Eq. (8) gives a poor result for most exeited vibrational levels of the 
a^E+ potential with an rRMS exeeeding 4% in the range of Vq = 2 to 27. This behavior 
is related to the faet that Franek-Condon faetors mainly depend on the amplitude of the 
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Figure 8. Analytical fits of the dynamical polarizabilities of the a^H^ molecule in Vq = 0, 
N = 0 (a), and of the molecule in vx = 0, N = 0 (b). Frequency ranges excluded from 
the fits are indicated as shaded areas. At the scale of the figure, the original curves and the 
fitted curves are indistinguishable. 


excited wave function around both the inner and the outer classical turning points of the initial 
vibrational level studied. The shape of the a^E+ potential is quite different compared to the 
relevant excited states. Consequently, the inner turning point region of a given vibrational 
level Va will induce couplings to levels of excited molecular states with energies strongly 
different from those related to couplings induced by the outer turning point region. This 
effect mainly occurs for the excited potential wells as they are deep. Instead, the depth 
of the potential is not sufficient to create such a variation. 

Thus we added one more transition term to the ansatz of Eq. (8) in order to account 
for both the inner and outer part of the a^E+-^E+ transitions in the model. This reduced 
significantly the rRMS of the effective polarizabilities of the levels (see Fig. 9). We 
want to emphasize that such an interpretation gives a reasonable physical picture for most 
levels. However, for some levels like the deeply bound ones or those close to the dissociation 
limits, the three-effective-transition model is somewhat artificial and the effective parameters 
should be taken only as numerical parameters needed to easily obtain the corresponding 
polarizability. 

9. Conclusion 

We have studied the dynamical polarizability of the ^^Rb 2 molecule in the rovibrational 
ground state of the a^E+ potential. Calculations of both, the real and imaginary part are 
provided and we measured Re{a((o)} at A = 1064.5 nm. Our experimental and theoretical 
findings show good agreement. From our computed value of Im{a((o)} at this wavelength, 
we expect trapping times of the molecules on the order of seconds for lattice depths around 
50Er, which was confirmed by our observations. We also have investigated theoretically the 
dynamical polarizability of the singlet ground state X^E+. These results are interesting for 
future STIRAP transfer of Rb 2 to the corresponding rovibronic ground state. Furthermore, we 
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Figure 9. Relative root mean square (rRMS) values of the analytical fits to the numerically 
calculated dynamical polarizabilities. Blue closed (red open) circles correspond to fit results 
obtained with the ”two-effective-transitions” ansatz for the vibrational levels Vq of a^X;+ (vx of 
X*X;+), whereas the black closed squares stem from fits including a third effective transition 
to represent the dynamical polarizabilities with respect to Va, which significantly increases the 
accuracy (see text). The rRMS values for vx > 40 do not exceed 0.1% (not shown here). 


have introduced a simple analytical expression to parametrize the dynamical polarizabilities 
for all levels of both, a^E+ and states. By fitting this expression to the numerical results, 
we have extracted effective parameters, which can be used to reproduce Re{a((o)} of a given 
vibrational state with high fidelity. The precise knowledge of the dynamical polarizability 
enables accurate control of optical dipole potentials and therefore is of importance for future 
experiments with deeply bound Rb 2 molecules. 
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Appendix 

A.l. Polarizability of the Rb^ core 

As the two ionic Rb+ cores are only weakly perturbing each other, we consider the molecular 
core polarizability as twice the atomic Rb^ polarizability. First, we calculate the atomic 
polarizabilities at imaginary frequencies aRb(/Ci)) [55], which includes only the contribution 
of the valence electron. We subtract them from the values of Ref. [23], where the resulting 
differences represent the contribution of the core electrons, and thus the Rb+ polarizability. 
Following Ref. [23] this estimate assumes that the influence of the valence electrons on the 
core polarizability is negligible. 
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Figure 10. Computed transition electric dipole moments (TEDMs) for the main transitions 
from the a^21+ and the X*X;+ states to the states correlated to the 5s + 5p dissociation limit (a), 
the 5s + 4ii dissociation limit (b), the 5s + 6s dissociation limit (c), and the 5s + 6p dissociation 
limit (d). The individual curves can be assigned to the corresponding transitions as indicated 
on top of the graph using the numbers i and j. At large distances the TEDMs converge towards 
dss-ni which is equal to the atomic TEDMs multiplied by v/2 (see text). 


We use an ansatz similar to the one of Eq. (8) to model the eore polarizability with 
two effeetive transitions. This yields the effeetive frequeneies CO/{Inc) = 165912em^^ and 
362918em ^, and the effeetive dipole moments 2.72799a.u. and 1.40119a.u. We note that 
these two transition frequeneies are on the order of magnitude of the main transition in Rb^ 
and Rb^+, respeetively. The eore polarizability is a small eontribution slowly varying from 
17.9 a.u. at vanishing frequeney up to 18.1 a.u. in the optieal domain relevant here. For the 
statie ease, we ean eompare the result of our ealeulations to the value of 18.2 a.u. reported in 
Ref. [47] and find good agreement. 

A.2. Transition electric dipole moments 

For the sake of eompleteness we show in figure 10 the transition eleetrie dipole moments 
ineluded in the ealeulations of the dynamieal polarizabilities as funetions of the intemuelear 
distanees. The eorresponding numerieal values are given in the Supplemental Material 
[20]. Some of the data with respeet to the lowest transitions have already been reported 
in Refs. [25, 26]. The largest TEDMs in the range of the PECs eoneern the first exeited 
E+ and IT potentials [see Fig. 10(a)]. At large distanees the moleeular exeited eleetronie 
wave funetions beeome elose to the form [ 05 j( 1 ) 0 „/( 2 ) ± 05 j( 2 ) 0 „/( 1 )]/a /2 and therefore the 
TEDMs eonverge towards d^s-nh whieh is equal to the atomie TEDMs multiplied by a / 2 . 
We find dss-sp = 4.23 a.u. and d^s^cp = 0.347 a.u., in exeellent agreement with the values 
extraeted from the NIST database [44] (4.226 a.u. and 0.3531 a.u.), obtained by averaging 
the TEDMs eorresponding to 5s — npi ^2 5s — np 2/2 for n = 5 and 6, respeetively. 
This eonfirms the good quality of the present representations of the atomie eleetronie wave 
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